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0.1.1 Introduction 



Q 

The discovery of non-periodic solids has motivated the construction of numerous examples of 
aperiodic tiling models, and led to the systematic theory of cut and project or model sets, as 
they are now called, see [[l],||,|3|] for details and further references. Most combinatorial ques- 
tions of crystallography, such as sublattice or shelling structures, have a natural analogue for 
aperiodic systems (we will explain the technical terms in more detail later on). However, the 
£N| ■ traditional methods of crystallography do not apply. Fortunately, the most important systems 

possess a high degree of symmetry, which manifests itself in an intimate relation to algebraic 
number theory [Eh. This relationship has been exploited successfully to tackle combinatorial 
questions for quasicrystals 

Interestingly, this approach also simplifies the treatment of crystals. A fairly common fea- 
ture is the use of Dirichlet series generating functions, which emerges from the observation 



that many counting functions, when properly normalised, can be expressed in terms of mul- 
tiplicative arithmetic functions, see JTo| ] for background material. This leads to a systematic 
and unified approach which will be summarised in this article. Besides explicit results, the 
generating functions also allow a precise calculation of asymptotic properties. 

There are various other combinatorial problems, such as the determination of coordination 
sequences and coronae, the orbit structure under inflation, or general complexity considera- 
tions. The results are usually encapsulated in terms of generating functions, compare Jll| ] 
for background material. In this expository article, we will concentrate on problems that are 
connected with Dirichlet series generating functions, and only summarize other developments. 

To be able to explain the concepts in a simple fashion, we shall mainly use two examples, 
namely the triangular lattice, written as r — Z[£ 3 ] with £ 3 = exp(27ri/3) = (—1 + i-\/3 )/2, 
the set (ring) of Eisenstein integers, and the vertex set of the twelvefold symmetric shield 
tiling, see Fig. [!} We concentrate on explicit results for these examples, and refer to original 
sources for a more general exposition and for details on the asymptotic behaviour. 



0.1.2 Counting general sublattices 

Let us first consider the triangular lattice r = Z[£ 3 ], or, in fact, any planar lattice, i.e., the 
integer span of two vectors in the plane that are linearly independent over K. Our first question 
is for the number £ 2 ( m ) of sublattices J" of given index m — [r : r'] . This is a multiplicative 



Figure 1: Finite patches of the triangular lattice and the twelvefold shield tiling. 



function, i.e., £ 2 {mn) — £ 2 (rn) l^in) for m, n coprime. Due to multiplicativity, it is best 
encoded in a Dirichlet series generating function, where one obtains [||, Eq. (3.9)] 

, „ £ (m) , , , s T — r 1 1 

^2* = E^V = c*c*-i = 111 — -i — — t 1 ) 

m s J - J -l — p s 1 — p L s 

m—l P 

_ 1iij.ij.iiiil2ilil5j.13i. _18_ I 12 I 28 I 

1 " r 2 s "r" 3 s 4 s "•" 5 s ' 6 s ' 7 s ' 8 s 9 s ' 10 s 11 s 12 s ' ■ ■ ■ 

Here, £(s) = X)m=i m s i s Riemann's zeta function, p in the product runs over all rational 
primes, by which, following the usual convention, we mean the primes of Z. Furthermore, 
£ 2 {jri) — Sdlm d * s tne ordinary divisor function, compare 0] for details. 

The corresponding question is well defined also for L = Z[£ 12 ],£ 12 = exp(27ri/12), when 
considered as a (free) Z-module of rank 4. The number of full rank submodules of index to, 
^ 4 (to), is again multiplicative and has the generating function ^ App. A] 

F 4 (s) = E^ = CWC(*-l)C(*-2)C(*-3) (2) 

m—l 

= 1 4- J! I 40 , 155 , 156 , 600 , 400 , 1395 , 1210 , 2340 , 1464 , 
2 s 3 s ' 4= " r 5 s ' 6 s ' 7= "r 8 s 9 s ' 10 s ' 11 s 

where £ 4 (m) = J2d ■ -d =m d i 4 ^4 with the sum running over positive di only. The 
answer is the same for all Z-modules of rank 4. They include the cyclotomic rings Z[£„] 
with n £ {5, 8, 12}, which are important for quasicrystals [|T2"|]. Note that we follow the 
mathematical convention to restrict to n ^ 2 (4) because Z[£„] = Z[£ 2j J for n odd. 

Clearly, this combinatorial question is purely algebraic, and the method also applies to 
counting finite index subgroups of finitely generated free Abelian groups, see App. A] for 
details. The generating function for counting index to subgroups of Z" reads 

F n( s ) = E^V 2 = C(*K(*-l)-...-C(*-n+l) (3) 

m—l 

and gives the formula £ n {m) = £ d . =m d i 4 ' ■ ■ ■ ' 
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0.1.3 Counting similarity sublattices 

Let us turn to (geo)metric properties and ask for the number a 6 (m) of triangular sublattices of 
r of a given index m, which coincides with the number of ideals of Z[£ 3 ] of norm m because 
Z[£ 3 ] is a principal ideal domain [Q. The answer is given by the Dedekind zeta function of 
the cyclotomic (also, quadratic) field Q(£ 3 ) = Q(V~ 3)» i- e -' one obtains [|l4], Eq. (9)] 



f" "\ V* u 6\" 1 ) 1 TT 1 TT 1 

m=l P=l(3) ^ ^ j PS2(3) 



= 14 — 4 — 4 — 4- — 4- — 4- — 4- — 4 — 4 — 4- — 4 

J- T 3 s '4 s 7 s ' 9 s ~ 12 3 ' 13 s ~ 16 s ~ 19 s ~ 21 s ~ 25 s ~ • • • 

where, e.g., p = 1(3) means that the corresponding product runs over all rational primes 
congruent to 1 modulo 3. 

The corresponding question for the module Z[£ 12 ] is answered by the Dedekind zeta func- 
tion ^ of the cyclotomic field Q(£ 12 ), which reads @ Eq. (13)] 

x - q 12 to 

m— 1 

4 s 9 s 13 s ~ 16 s ~ 25 s ' 36 s ' 37 s ' 49 s ' 52 s ~ 61 s ' " " ' 

These two generating functions also have an interpretation in terms of planar colourings 
|[l4|,[jl], p7j]: a 6 (m) is (up to permutation) the number of colourings of the triangular lattice 
r where one colour occupies a similarity (hence triangular) sublattice of r of index to and 
the remaining colours occupy its cosets. The function a 12 (m) then counts the colourings of 
Z[£ 12 ] via similarity submodules of index m, all of which are principal ideals of Z[£ 12 ] as the 
latter is a principal ideal domain, see [p3[p|] for details. For a discussion of the corresponding 



colour groups, see 1 16] 



This approach can be extended to all cyclotomic fields with class number one (the corre- 
sponding rings of integers are then principal ideal domains), see []l5|,[l4|^8|], and references 
given there, for details. 

0.1.4 Counting coincidence sublattices 

Another geometric problem, with interesting applications to grain boundaries and twinning 
phenomena in crystals and quasicrystals, is the classification of sublattices of r that can be 
seen as the intersection of r with a rotated copy of itself. These are the coincidence sublat- 
tices, see Q| for details, and the corresponding coincidence index is called Z'-factor in mate- 
rials science. Once more, the generating function (for the number of coincidence sublattices 
of index to) is best written as a Dirichlet series, where, using r — Z[£ 3 ], one finds 

— I J 2 I 2 i 2 | 2 | 2 | 2 | 2 | 2 , 2 , 2 , 

7 s 13 s ' 19 3 ~ 31 s ~ 37 s T 43 s ' 49 s ~ 61 s ~ 67 s ' 73 s ' • • ■ 



which is derived in [g, Sec. IV. A]. 

Behind this combinatorial problem is a group structure. It turns out that the set of rotations 
R such that r n RT is a coincidence sublattice of r forms a group, called SOC(-T), which 
is isomorphic to C e <S> Z^*°\ see [ |l9[ ] for details, in particular on the structure of this Abelian 
group. Also, the possible coincidence indices [r : (r n Rr)] =: S(R) (excluding and oo) 
form a monoid (a semigroup with unit), generated by the (rational) primes p = 1 (3). 

The formulation with cyclotomic integers admits an extension of the previous results to 
Z[£ 12 ], where the generating function for the number of coincidence submodules of Z[£ 12 ] of 
index to is Sec. IV.F] 



$ (s) - n ( 1+P ~' V n 1+p " 2 ' - Cq(Ci2)(s) (?) 

pE=l(12) V ^ 7 p=±5(12) ^ S Q(\/3) V ' 

T 13' T 25 s ' 37 s ' 49 s ~ 61 s ~ 73 s ~ 97 s T 109 s ^ 157 s ~ ' " 



Here, Cn^yg^ 5 ) i s me Dedekind zeta function of the quadratic field Q(v3 ), i.e., 

^Q(V3)( S ) = l _ 2-s 1 _ 3-s II n _ p - S )2 II 1 _ „-2s ■ (8) 

p=±l(12) V y ' p=±5(12) 

There is one subtlety in the application of this result to a discrete point set such as the 
vertex set of the shield tiling. A small acceptance correction factor is needed, which can be 
calculated explicitly, except for situations where the window has fractal boundary, as in |[20||. 
Further details, together with a general discussion of the case of 71-fold symmetry, can be 
found in [Si. 



0.1.5 Central shelling 

Let us discuss the shelling problem, first in its version for the central shelling. Here, one asks 
for the number of points of T = Z[£ 3 ] on circles of radius r around the origin. The result is 
usually given in terms of lattice theta functions, see Ch. 4] for an extensive exposition. 
However, in our situation, it can also be encapsulated in a Dirichlet series. To this end, one 
considers only radii r > and divides the corresponding shelling number, c(r 2 ), by 6, which 
is the trivial symmetry factor. If a shell is non-empty, we have r 2 = x x with x £ r, and 
r 2 = rn is a (rational) integer. What then remains is the multiplicative function a 6 (to) whose 
generating function was given above, in a different context, in Eq. (^). Further details and 
references are given in [|2|] . 

In the case of Z[£ 12 ], the central shelling function c depends on r 2 £ Z[£ 12 + £ 12 ] = 
Z[V3], and one needs the primes of Z[%/3 ] (see Jl3| ] for details) to derive a formula for c(r 2 ), 
compare [Jjl^]. Let us write p for a prime in Z[v3], and let t(p) be the highest power t such 
that p *■ divides r 2 , written as p 1 1| r 2 . Then, whenever t{p) is odd for a p that is also a prime in 
Z[£ 12 ] (i.e., an inert prime), one has c(r 2 ) = 0. Otherwise, one finds the formula 



c(r 2 ) = 12 J] (*(P) + 1) 



(9) 
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z 


Z[V3] 


Z[fu] 


2 


2 = (2-V3)-(l + V3) 2 


1 + V3 prime 


3 


3=(V3) 2 


Vo prime 


p=l (12) 


p = pp' 


p = PP, p' = P'P 7 


p = -1 (12) 


p = pp' 


p,p' prime 


p = ±5 (12) 


p prime 


p = PP 



Table 1: Splitting structure of primes for Z[£ 12 ]. Here, ' is a Galois automorphism of Q(£i2) 
which acts as v3 i— > — \/3 in Q(V3 ), and denotes complex conjugation. Note that 2 — v3 = 
(2 + \/3 ) _1 is a unit in Z[\/3]. 



where the product runs only over a representative set of those primes p of Z[\/3] that split 
in the extension to Z[£ 12 ]. In this example, they are precisely the primes that originate from 
rational primes p = 1 (12) and p = ±5 (12) . 

Fortunately, the function c(r 2 ) need not be calculated for all possible values of r 2 . Since 
the units of Z[\/3] are of the form ±(2 + \/3) with I S Z, all positive units are actually 
totally positive, i.e., also their algebraic conjugate in Q(-\/3) is positive. The fundamental 
unit is 2 + V3 = (1 + £12X1 + £12)' nence a H positive units are of the form u — xx with 
x e Z[£ 12 ], and this implies that c(r 2 ) only depends on the (principal) ideal of Z[\/3] that 
is generated by r 2 , compare the proof of theorem 1 in [^J. So, whenever s 2 = ur 2 with u 
a positive unit, we have c(s 2 ) = c(r 2 ). Also, one can see that c(r 2 ) = c((r 2 )'), where ' 
denotes algebraic conjugation in Q(\/3 ), defined by the Galois automorphism \/3 — \/3, 
compare Table []]. 

Can the result be given in an even simpler way? Suppose, for a moment, that the function 
c(r 2 ) would only depend on the norm of r 2 , i.e., on m = r 2 (r 2 )', where to is a rational in- 
teger. This would mean that c(r 2 ) — 12/(to) with /(to) a multiplicative arithmetic function, 
whose Dirichlet series turns out to be 

/M 



CM = do) 



m— 1 

1 1 



1 - 4-s 1-9 



^1(12) 



n ^— n — 

11 l_p-2 S 11 d_ D -2 S \ 
p=-l(12) ^ p=±5(12) ^ P I 



= 1 4- -L 4- -L J ?^ I 1 I I 1 I 1 I 1 I 1 I 1 L 

4 s 9 s 13 s 16 s ~ 25 s T 36 s T 37 s ' 49 s T 52 s T 61 s ' ' ' ' 

As can be seen from the exact shelling formula (^), this cannot provide the general answer. 
However, it is correct for many cases, and is then by far the simplest way to calculate c(r 2 ). 

The first non-trivial value of m where this fails is to = 169 = 13 2 . Observing that 
13 = (4 + \/3)(4- V3) = pp' in Z[v3], the possible values of r 2 with norm m are p 2 , 
pp' and (p') 2 , with shelling numbers 36, 48 and 36, respectively. Since /(169) = 3, we get a 
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Figure 2: Vertex set of the shield tiling and its dodecagonal window (black). Also shown are 
three shifted copies of the window (grey), see text for details, as needed for Table ^. 

wrong answer for the middle case this way. This is an example with a prime p = 1 (12). The 
contributions to non-empty shells from the other rational primes are properly counted by the 
Dirichlet series (|l0|) because they either correspond to a unique prime in Z[\/3] or to primes 
which do not split in the step from Z[V3] to Z[£ 12 ], see Table [l| 

So, let us assume that r is the radius of a non-empty shell, i.e., t(p) is even for all inert 
primes p of Z[\/3]- Then, 12/(m) with m — r 2 (r 2 )' gives the correct shelling number if, 
whenever m is divisible by a rational prime p = 1 (12), r 2 is divisible by either p or p', 
compare Table |], but never by both. Beyond this situation, one has to use formula (^). 

Eq. (^), together with the partial simplification of Eq. (|To| ), gives the result for the central 
shelling of the full module. When passing to a model set ||l|], which is a discrete subset of the 
full module, an obvious selection takes place for the possible radii. Furthermore, a correction 
factor may become necessary which depends on the window chosen, see Jl8| ] for details. A 
systematic approach to a complete generating function would employ summatory functions of 
arithmetic characters [p5j], which is postponed to future work. 



0.1.6 Averaged shelling 

In this last application, we are interested in the average number of points on shells of radius 
r, where the average is taken over all points of our point set as possible centres. In the lattice 
situation, the averaged and central shelling numbers coincide, but this is not so for model 
sets [E2JII]. To be specific, we explain this for the vertex set of the shield tiling. As for model 
sets in general, all averages needed are well defined and unique. 

Let us first describe the vertex set of the shield tiling [ ^3| ] in algebraic terms. We use L — 
Z[£ 12 ] as above, together with the Galois automorphism * of Q(£ 12 ) given by £ 12 1— > (£i 2 ) 5 - 
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r 2 


representative 


orbit length 


norm 


shift type 


a(r 2 ) 


2 


-73 




12 


1 


2 


8-2\/3 


4- 


- 2V3 


€ + e - 1 


12 


4 


1 


2 


6- 


- 3a/3 




12 


9 


2 


4-2V3 




1 


i 


12 


1 


1 


8 


5- 


- 2^3 




24 


13 


3 


10-4V3 




2 


1 + e 3 


12 


4 


2 


48 - 24 V3 


4 


-V3 


e + e 2 - 2 


24 


13 


3 


6 


8- 


- 3V3 


3e - e - 1 


24 


37 


3 


-76 + 44\/3 




3 




12 


9 


1 


-4+f V3 


7- 


- 2V3 


2 - 2£ - £ 2 


24 


37 


3 


20- f ^3 


2 


+ V3 


1 + C 


12 


1 


2 


48 - 22v^ 




4 


2 


12 


16 


1 


2 + 1^3 



Table 2: Averaged shelling numbers a(r 2 ) for the shield tiling of edge length 2sin(7r/12) = 
V 2 — for all possible distances < r < 2. Representative elements (with f = f 12 ) 
and orbit lengths refer to the point symmetry group D 12 of the tiling. The norm of r 2 is m = 
r 2 (r 2 )', and the shift types are explained in Fig.0. 



This is the star map of the standard model set construction [|l|], which gives 

A = {x G L | x* G W} (11) 

where the window W is a relatively compact set with non-empty interior. For the shield tiling, 
we choose W as a regular dodecagon of edge length 1, hence of inradius R K = (2 + \/3 )/2 
and circumradius i? > see Fig. ^]for the correct orientation. 

If the window is centred at 0, the model set is singular, while generic examples are obtained 
by shifting the window, i.e., A u = {x E L | x* G W + u}, for almost all u £ M 2 . The 
shortest distance between points in a generic A u is (vo — l)/v2. Joining such points by 
edges results in a shield tiling; different generic choices of u result in locally indistinguishable 
shield tilings. 

Table [| lists all distances < r < 2 between points of A which occur with positive fre- 
quency, together with a representing point in L. The corresponding orbit length is given w.r.t. 
the point symmetry group D 12 . Here, it coincides with the value of the central shelling number 
c(r 2 ) — 12/(m), as calculated by means of Eq. (||) or Eq. (|T(i|). For larger values of r 2 , c(r 2 ) 
can comprise several Z? 12 -orbits, but only complete ones because the window is twelvefold 
and reflection symmetric. For the averaged shelling, we need to calculate the relative overlap 
area of the window W and a shifted copy (the so-called covariogram of W), where the shift is 
the star image of the representing vector Jl8|]. For our examples, all three relative orientations 
indicated in Fig. || emerge, with squared midpoint distance (r 2 )' in Q(-\/3 )• The resulting 
averaged shelling number (last column) is always an element of Q(y3). More precisely, all 
averaged shelling numbers are in |Z[-\/3]. This follows from the fact that the frequency mod- 
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ule of the vertex set of the shield tiling is ggZ[\/3]- The frequency module is the Z-span 
of the frequencies of all possible finite patches and can be determined from the topological 
invariants of the shield tiling . 

For the overlap scenarios 1 and 2 of Fig. |[ where the shifts are along the two principal 
reflection axes of the dodecagon, the covariogram has a relatively simple form. We find 



his) 



15+2V3 

18 

5+2V3 _ 
6 

0, 



2V3 . 



2V3-3 2 
18 * ' 



2V3-3 2 
6 * ' 



< s < 1 

1 < s < 2R K - 1 
2-R< - 1 < s < 2R< 
2R K < s 



(12) 



and 



h 2 (s) = < 



V2(x/3-l) 



5+V3 _ 
6 

4+2y/3 
3 



3 6 1 6 
V2(l+V3) 



7-4V3 2 
6 * ' 



is 2 



< s < R > 
R > <s<R > +V2 
R>+V2<t< 2R > 
2i? > < s 



(13) 



These functions have to be multiplied by the orbit length to give the averaged shelling number, 
resp. the corresponding contribution to it. 

In Table ^J, the four examples of type 3, where the shift is not along a symmetry direction, 
were calculated separately (and exactly). We do not give the more involved general formula 
for the covariogram of the regular dodecagon, but mention that 



h(s) = — arccos , . 

^ ' 7T \2R) 



ttR 



i-f-V 

\2RJ 



(14) 



with R = R > , is a very good approximation. It is obtained by replacing the dodecagon 

W with a disk of equal area, hence of radius R. This function is also known as Euclid's hat, 
see J25| , p. 100] and references given there for details. 

The averaged shelling in one dimension is rather simple and can be given in closed form 
"J], though no appropriate approach via generating function has been formulated so far. In 
two dimensions, the example of the Ammann-Beenker rhombus tiling with eightfold symme- 
try is treated in [ p2"[|T8[ |, while examples with tenfold symmetry are shown in [^2j| (rhombic 
Penrose tiling) and pq ] (Tubingen triangle tiling). So, with the example shown here, the stan- 
dard symmetries in the plane are covered. Again, the understanding in terms of generating 
functions is missing, and in three or more dimensions, the averaged shelling has not been 
looked at thoroughly so far. 



0.1.7 Other combinatorial questions 

Another combinatorial problem, which can be seen as a variant of the shelling problem, is 
the determination of the coordination structure [EU271E8]. This needs a graph structure in 
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addition to the point set, and is usually formulated for tilings. Once again, one can ask for 
the number of points in the nth coordination shell (or corona), or for the analogous averaged 
quantity. The latter is more suitable for non-periodic (face to face) tilings. If they are con- 
structed by a primitive inflation rule, or by the standard cut and project method (e.g., in its 
dualization version [^9[|), we know that all averages exist and the problem is well posed. 

Let us first consider the situation of a lattice. If we denote the number of points at graph 
distance n by s n , one uses the ordinary power series generating function S(x) = X^^Lo s n a; ™- 
For the class of root lattices and root graphs, these generating functions are rather simple and 
can be given in closed form, see [[30, 27]. For aperiodic tilings, one faces the same difficulties 



as for the averaged shelling discussed above. The rhombic Penrose tiling and its octagonal 
sibling, the Ammann-Beenker tiling, are discussed in [pq], Once again, the averaged coordi- 
nation numbers seem to be "nice" numbers. This is due to the structure of the corresponding 
frequency modules, and can be quantified [g4J]. 

Closely related is the more general question for the patch counting function, which is 
a direct measure of the local complexity of the structure. For example, it is well known 
that Sturmian sequences possess n + 1 different words of length n. This is the smallest 
complexity possible for any non-periodic structure in one dimension. To what extent this can 
be generalized to higher dimensions is still unsettled, see for details. In order to come to 
a natural analogue to the statement about the Sturmian sequences, one could make use of the 
knowledge obtained from the shelling function discussed above. It provides the possible shell 
radii and hence the sizes of the patches to be considered. 

The method of iterated inflation is, besides the projection method, the most frequently used 
approach to generate aperiodic tilings with long-range order. Under some mild conditions, the 
resulting structure is linearly repetitive, see [[$2]] for details on this concept. Inflation also 
acts naturally on the set of all tilings that are locally indistinguishable (LI) to the fixed point 
constructed, i.e., on the entire LI class. In view of the high degree of (local) repetitivity, one 
is also interested in the orbit structure of the LI class under the inflation /. Let a n denote the 
number of fixed points of I n , n > 1, and define 



CuCO = exp V^z" 1 (15) 




so that Z 5 LI /^ = Yl^-i a m z m is the usual power series generating function of the sequence 
{a m | m > 1}. The object £li( z ) i s called a dynamical zeta function, see [ p3| ] and references 
given therein for background material. 

This function can be calculated systematically, as explained in [Q. The cycle structure 
of inflation can also be extracted, e.g., by means of the Euler product expansion 

where c n is the number of n-cycles under inflation. Clearly, a n = X) m |n mCm ' an( ^ an 
inversion is possible employing the Mobius function, compare JlO|]. 

If the tiling admits a construction by the projection method, it also admits a universal torus 
parametrization |B5lR6||. The latter will be one-to-one on generic members of the LI class, but 
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Fibonacci 


Penrose 


LI class 


1 - z 


(l-z - z 2 ) 2 (l + z) 


1 — z — z 2 


(l-3z + z 2 )(l + z-z 2 ) 3 (l-z) 


torus 


l-z 2 


(l-z-z 2 ) 2 (l + z-z 2 ) 2 


1 — z — z 2 


(l-3z + z 2 )(l-z) 2 (l + z) 4 



Table 3: Dynamical zeta functions for the inflation action on the Fibonacci chain and the rhom- 
bic Penrose tiling, both on the full LI classes and on their torus parametrizations. 



multiple-to-one on singular members. Consequently, since inflation is compatible with this 
reduction step, there is another dynamical zeta function, this time for the orbit structure on the 
torus. Both types of zeta functions are known to be rational. Two examples are summarized 
in Table 9, further examples can be found in B4l||,|l. 



0.1.8 Open problems 

The above examples demonstrate that a systematic and unified approach to combinatorial 
problems of (quasi)crystallography is possible, at least in the planar case. The situation is 
more involved in dimensions d > 3, where satisfactory results so far exist only for the most 
symmetric cases, i.e., those with (hyper)cubic or (hyper)icosahedral symmetry [RH,5^,^,§]. 
Still, even some of these cases leave room for improvement and simplification, e.g., along the 
lines mentioned around Eq. (|Io|). 

On the other hand, model sets with high symmetry are closely related to lattices in higher 
dimensions, where many of these questions are still open, compare plj]. One may expect 
some progress at least for the class of root lattices, hence also for quasicrystals derived from 
the m Q39 ]. Further progress is also needed in the investigation of colour symmetry groups, 
see Jl6|] for a summary of the present state of affairs. 

A big mystery is the meaning of the averaged shelling function. For the standard tilings 
with all magic properties (inflation rule, perfect matching rules, pure point diffraction etc.), the 
averaged shelling numbers always seem to be "nice" (being algebraic integers or rationals with 
bounded denominator), while this is not the case for model sets with generic windows. This 
phenomenon points towards another function defined by these numbers, with various analytic 
and topological consequences, but we do not know how to substantiate this at present. 
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